## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

## Warning in geom_rect(aes(x = TYPE, xmin = as.numeric(as.factor(TYPE)) - :
## Ignoring unknown aesthetics: x

## Linking to GEOS 3.11.0, GDAL 3.5.3, PROJ 9.1.0; sf_use_s2() is TRUE
## Loading required package: viridisLite
## Classes 'sf' and 'data.frame':   170 obs. of  4 variables:
##  $ MODZCTA           : num  10001 10002 10003 10026 10004 ...
##  $ median_price      : num  570 831 893 595 1036 ...
##  $ num_establishments: int  70 192 77 49 10 10 7 21 76 46 ...
##  $ the_geom          :sfc_MULTIPOLYGON of length 170; first list element: List of 1
##   ..$ :List of 1
##   .. ..$ : num [1:84, 1:2] -74 -74 -74 -74 -74 ...
##   ..- attr(*, "class")= chr [1:3] "XY" "MULTIPOLYGON" "sfg"
##  - attr(*, "sf_column")= chr "the_geom"
##  - attr(*, "agr")= Factor w/ 3 levels "constant","aggregate",..: NA NA NA
##   ..- attr(*, "names")= chr [1:3] "MODZCTA" "median_price" "num_establishments"
##  - attr(*, "na.action")= 'omit' Named int [1:4] 83 135 137 147
##   ..- attr(*, "names")= chr [1:4] "83" "135" "137" "147"
## Rows: 2639 Columns: 17
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (10): County, Operation Type, Establishment Type, Entity Name, DBA Name,...
## dbl  (5): License Number, Zip Code, Square Footage, Counties, New York Zip C...
## lgl  (2): Address Line 2, Address Line 3
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
## Warning in inner_join(retail_food_stores, ny_house_data, by = c(`Zip Code` = "Main_address_zip")): Detected an unexpected many-to-many relationship between `x` and `y`.
## ℹ Row 1 of `x` matches multiple rows in `y`.
## ℹ Row 80 of `y` matches multiple rows in `x`.
## ℹ If a many-to-many relationship is expected, set `relationship =
##   "many-to-many"` to silence this warning.
## `summarise()` has grouped output by 'Zip Code'. You can override using the
## `.groups` argument.

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ forcats   1.0.0     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
## 
## Attaching package: 'plotly'
## 
## 
## The following object is masked from 'package:ggplot2':
## 
##     last_plot
## 
## 
## The following object is masked from 'package:stats':
## 
##     filter
## 
## 
## The following object is masked from 'package:graphics':
## 
##     layout

Load dataset